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ABSTRACT 

We report measurements of the Diffuse Galactic Light (DGL) spectrum in 
the near-infrared, spanning the wavelength range 0.95-1.65 /mi by the Cosmic In¬ 
frared Background ExpeRiment (CIBER). Using the low-resolution spectrometer 
(LRS) calibrated for absolute spectro-photometry, we acquired long-slit spectral 
images of the total diffuse sky brightness towards four high-latitude fields spread 
over four sounding rocket flights. To separate the DGL spectrum from the total 
sky brightness, we correlated the spectral images with a 100 /mi intensity map, 
which traces the dust column density in optically thin regions. The measured 
DGL spectrum shows no resolved features and is consistent with other DGL mea¬ 
surements in the optical and at near-infrared wavelengths longer than 1.8 /mi. 

Our result implies that the continuum is consistently reproduced by models of 
scattered starlight in the Rayleigh scattering regime with a few large grains. 
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Subject headings: infrared: ISM — ISM: dust, extinction — ISM: general - 
scattering 


1. Introduction 


Diffuse Galactic Light (DGL) arises from stellar radiation scattered by dust in the 
interstellar medium. The DGL spectrum includes information on the optical properties of 
interstellar dust, such as the grain size distribution and composition, as well as the interstellar 
radiation field (ISRF). 


Historically, DGL was originally d etected at optical wavelengths and interpreted as 


starlight scattered by interste l lar dust (lElvev &; Roach! 119371: iHenvev fe Greenstein 


van de Hulst fe de Jond 119691: lElsasser fe Haud 11960 ; 


Wolstencroft fc Rose 


1966; 


1941 


Mattila 


19791 1 . The interstellar dust has been studied through its emission properties in the far- 


infrare d fide Vries fe Le Poolel Il985l: lLaureiis et al.l 119871: IStarkl 119951 1 , and was mapped by 
IRAS ( Low et al. 1984 L 


In the optically thin limit, DGL scales with the light absorbed and re-radiated in the 
far-infrared, with a proportionality that depends on the properties of the dust grains. This 
linear correlation between DGL and far-infrared intensity enables us to separate the DGL 
component from other diffuse emission, such as Zodiacal Light (ZL), sun light scattered 
by interplanetary dust, Integrated Star Light (ISL) from undetected stars, and the Extra- 
galactic Background Light (EBL). 


A wavelengths shorter than 0.8 /mi, DGL has a significantly shallower spectral slope th an 
the ISRF due to the wavelength-dependent scattering cross-section of the grains (Brandt & Draine 


2012). This suggests that the DGL spectrum depends on the size distribution of the grains. 


The scattering cross-section can be approximated by Rayleigh scattering theory, where the 
size of the dust grains is smaller than the wavelength. Near-infrared measurements may be 
more sensitive to the size distribution of dust grains than measurements at optical wave¬ 
lengths. 


Although interstellar dust scattering has been studied in several measurements, the size 
distribution is still under discussion. Kim & Martin (1199511 insisted that the size distribution 
requires grain sizes ranging from 0.003 /mi to 3 /mi with a peak at 0.2 /mi based on the 
polarization of star light assuming spheroidal dust particles. Weingartner & Draine! (120011 ). 
hereafter WD01, compared their own model with the observed extinction of starlight, and 
claimed that interstellar dust includes l arge grains with radii a > 0.2 /mi with a half-mass 
radius ao .5 = 0.12 /mi where Draine (2011) defines ao .5 as 50% of the mass in grains with 
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a > ao. 5 - A model of Zubko et al. 020 041 ). hereafter ZDA04, consists of a small population 
of large grains with radii a > 0.2 /mi and many small grains, giving a half-mass radius 
a 0 .5 = 0.06 /im. 

Measurement of the DGL spectrum in the near-infrared helps determine the size dis¬ 
tribution of interstellar dust. However, DGL has not been measured from the ground at 
these wavelengths because DGL is low-surface brightness and much fainter than atmo¬ 
spheric airglow emission that contaminates large spatial scales. Even at optical wavelengths, 
ground-based DGL measurements are problematic, and suffer from systematic error due to 
airglow emission. Thus we measure the near-infrared DGL from sp ace us ing the sou nding 
rocket-borne Cosmic Infrared Background ExpeRiment (GIBER) (Bock et ah 20060 . We 
correlate measured brightness in the near-infrared with diffuse thermal dust emission in the 
far-infrared to extract the DGL measurement. Our results are the first DGL measurements 
at these wavelengths. 


2. Cosmic Infrared Background Experiment 
2.1. Low Resolution Spectrometer 

CIBE R is designed to stu dy diffuse near-infrared emission from abo ve the Earth’s at- 


mosphere ( 

Zcmcov et al. 

2013) 

Korneut et al. 

2013), including 


the spectrum of diffuse light in 0.8< A < 1.8 /mi ( Tsumura et al. 2013c). The LRS consists 


of an optical collimator which brings an image of the sky to focus on a mask containing 5 
slits, each spanning a field of view (FOV) of 5°x2.7 / . The light passing through the slits is 
then dispersed by a prism and brought to focus again on a 256x256 pixel HgCdTe detector 
array. 

The LRS covers the wavelength range 0.8-1.8 /mi with a resolving power R = AA/A = 
15-30. The FOV provides a large etendue to measure the diffuse light with high signal to 
noise ratio (S/N). To measure the dark current of the detector, a cold shutter cooled to 77K 
is mounted just before the detector, and is closed before and during the flight. In addition, 
the slit mask provides a masked region where light does not fall on the detector, to monitor 
any short-term changes in dark current. The 3a sensitivity of the LRS with 200 spatial pixels 
in a 50 sec integration is 2.5 nW m~ 2 sr _1 at 1.25 /mi, which enabl es ac curate measurement 
of the DGL brightness. Details of the LRS design are described in Tsumura et ah (12013d )). 



























Table 1. 


Our observed fields of the second flight and the fourth flight. RA and DEC indicate equatorial longitude 
and latitude, respectively, while 1 and b indicate galactic longitude and latitude, respectively. 


Field Name 

Exposure time [sec] 

Altitude [km] 

(RA, DEC) [degree] 

(1, b) [degree] 

Second flight 

SWIRE ELAIS-N1 

83 

139-172 

(243.069, 55.283) 

(84.89, 44.62) 

North Ecliptic Pole (NEP) 

67 

199-220 

(270.871, 66.004 

(96.13, 29.81) 

ElatlO (/? = 10 degree) 

9 

245-310 

(234.337, -8.466) 

(356.88, 46.08) 

Elat30 (/? = 30 degree) 

18 

314-295 

(223.058, 20.658) 

(23.52, 63.31) 

BOOTES-A 

63 

288-232 

(218.806, 35.120) 

(58.76, 66.79) 

Fourth flight 

DGL field 

65 

272-401 

(251.97, 68.85) 

(100.37, 36.17) 

North Ecliptic Pole (NEP) 

60 

425-505 

(270.82, 66.24) 

(96.02, 29.48) 

Lockman Hole 

50 

520-558 

(161.23, 58.58) 

(149.07, 51.65) 

ElatlO (/3 = 10 degree) 

50 

566-557 

(190.49, 8.02) 

(295.80, 70.77) 

Elat30 (/3 = 30 degree) 

50 

577-562 

(193.05, 27.96) 

(111.34, 89.15) 

BOOTES-B 

55 

555-509 

(217.23, 33.18) 

(54.90, 68.13) 

SWIRE ELAIS-N1 

55 

395-275 

(242.84, 54.77) 

(84.56, 44.64) 




Note. 


The coordinate systems are based on J2000. 
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Fig. 1.— Example of an LRS spectral image towards the north ecliptic pole. The left panel 
shows the processed flight image used to derive the DGL spectrum. The five vertical sections 
correspond to the locations of the spectrometer’s slits dispersed in the x direction. The y 
direction has spatial information. Bright stars detected as point sources appear as discrete 
stripes. The top, bottom and left edge of the image are masked, to monitor any short-term 
changes in dark current. The right panel shows the image after subtracting dark current and 
masking bright stars to isolate the diffuse emission. 
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2.2. Observations 

CIBER was flown four times, in February 2009, July 2010, March 2012 and June 2013. 
The payload was successfully recovered and refurbished after the first three flights. For 
the fourth flight, a larger launch vehicle was implemented which resulted in a much higher 
altitude apogee, but this configuration did not allow the payload to be recovered. The 
observed fields are listed in Table [T| with exposure time and altitude. The first three flights 
used a two-stage rocket launched from White Sands Missile Range in New Mexico, USA. The 
apogee on these flights was typically 330km, providing a total exposure time of ~ 240 seconds. 
In the fourth flight, the rocket was launched from Wallops Flight Facility in Virginia, USA, 
using a four-stage rocket. The payload reached 550 km with a total exposure time of 335 
seconds. For this study we use only data from the second and fourth flights, as the first 
flight data were contaminated by excess stray thermal radiation from the rocket skin, and 
in the third flight the LRS was operated as a polarimeter. 

The raw data, which are non-destructively sampled by the integrating detectors, were 
telemetered to the ground from the rocket during the flight. The celestial attitude control 
system achieved a pointing stability of < 8". 


2.3. Field Selection 


For the DGL analysis, fields with large contrast in iioo/mi emission across an LRS FOV 
are selected. These include the NEP field, observed in both the second and fourth flights, a 
field referred to as DGL, specifically targeted for its large expected dynamic range in DGL, 
and the ElatlO field, o bser ved in the second flight. Figures [2] and [3] show the 100 //rn intensity 
map flSchlegel et al.l 19981 ) of all the fields measured in the second and fourth flights. The 
NEP and DGL fields are high ecliptic latitude fields (see fourth column of Table [T]) and have 
low ZL brightness. We also select low ecliptic latitude fields with bright ZL, ElatlO and 
Elat30, in order to make a template spectrum for ZL separation in the observed data. 

In all the flights, there is emission from terrestrial atmospheric airglow, exospheric 
atmospheric airglow, and dissociated water vapor that outgasses from the payload early in 
the flight. Although airglow emission dominates the observed spectrum from 1.5 /im to 
1.7 /im in the first field of every flight, airglow emission decays exponentially with time 
and altitude flTsumura et all 2Q10) and is negligible for the remaining fields. Since the 
airglow emission is assumed to have no spatial fluctuation in the FOV of the LRS, we do 
not subtract it in our DGL measurement. We discuss the systematic error due to airglow 
emission in Section [0 
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Fig. 2. — The color maps indicate the 100 /iin intensity of the far-infrared cirrus emission 
flSchlegel et ah 1998 ) for the second flight. The gray lines present the field of view of the 
LRS. The equatorial coordinate systems are based on J2000. 
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Fig. 3.— As Figure [2] for the fourth flight 
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3. Instrument Calibration 
3.1. Wavelength Calibration 



Fig. 4.— Relation between wavelength of the incident light and the position of the detector 
pixels for the second flight. 


In order to measure the electromagnetic spectrum of diffuse light, it is essential to cal¬ 
ibrate the wavelength response, which encodes the mapping of each wavelength of incident 
light to a position on the detector array. This is determined through a series of laboratory 
measurements. Two different light sources are used for spectral calibration consisting of the 
SIRCUS _{the Spectral Irradiance and Radiance Calibrations using Uniform Sources facility, 


Bro wn et ah (120061 )) laser facility and a standard quartz-tungsten-halogen lamp coupled to 


a monochrometer. The wavelength of the SIRCUS laser is determined with an external 
wavemeter calibrated by National Institute of Standards and Technology (NIST). The wave¬ 
length of the monochrometer is calibrated using a He-Ne laser and the spectral lines of a 
Ne lamp. In both cases, monochromatic light was coupled to an integrating sphere through 
a fiber which illuminated the LRS aperture. Following each exposure with a monochromic 
light source, we fit the detected signal with a Gaussian function. The center of this Gaus¬ 
sian, in combination with the externally determined wavelength of incident light, is used to 
generate the wavelength map shown in Figure [H The wavelength accuracy of the input light 
is better than lnm, sufficiently small compared to the > 10 nm wavelength resolution of the 
LRS. 
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wavelength [/itn] 


Fig. 5.— Comparison of the calibration factor obtained with the SIRCUS laser and the 
super-continuum laser. The red line indicates the calibration factor of super continuum 
laser. The black dots give the calibration factor of the SIRCUS laser. 


3.2. Surface Brightness Calibration 


Absolute spectro-photome tric calibration is done through laboratory measurements, us¬ 
ing the techniques described in Zerncov et ah] (20131). Two different light sources are utilized 
in this calibration measurement: a super-continuum laser (SCL) for broad-band measure¬ 
ment, and the SIRCUS laser facility. These light sources are coupled to a 4” diameter 
aperture integrating sphere 0 whose port is viewed by the LRS. The absolute brightness 
of the integrating spheres viewed by the LRS is determined using absolutely calibrated ra¬ 
diometers and a monitor detector. During the measurements, the light sources are shuttered 
and the measured ambient signal is subtracted from the data. 


The dynamic range mismatch between the radiometer and the LRS requires a two-step 
bootstrapping approach. To avoid non-linearity effects in the LRS detector array introduced 
at high-photo current levels, the brightness of the integrating sphere must be attenuated. 
The absolutely calibrated radiometers measure the absolute brightness of the integrating 
sphere with light levels 10 4 ~ 10 6 brighter than for the LRS measurement. The monitor 
detector measures the intensity of a secondary smaller injection sphere with sufficient S/N 
in both cases to measure the coupling factor between the two measurement points. Therefore, 
when the intensity is reduced to levels accessible to the LRS, the absolute brightness of the 
main sphere is inferred using the monitor signal and the previously determined ratio. The 


1 Manufactured by Gigahertz-Optik, Inc. 
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calibration uncertainty of the radiometers is approximately 0.3 %, as quoted by NIST. 

We calculate the calibration factor (CF), ^ Ix ^ x ~ x from both the SIRCUS laser and 
SCL measurements, where I\ indicates the absolutes brightness of the integrating sphere 
viewed by the LRS, /(A — A') is the response function of the LRS d etermined from the slit 
width and the point spread function measured in the laboratory (Tsumura et a l. 2013d). 
and i is the photocurrent measured by the LRS during the calibration measurement after 
correcting for the non-linearity of the detector flTsnmnra et al. 2013d ). 


Figure [5] shows the CFs for the second flight after correcting for the transmittance of 
the additional window used in the laboratory. The la statical uncertainty is estimated to 
be < 0.1 % from the variance across all of the detector array pixels. The measured CFs are 
consistent within a 3 % r.m.s variation, which sets our systematic uncertainty of the surface 
brightness calibration. 


3.3. Flat Field Correction 



x-pixel 


Fig. 6.— Flat-field correction map of the second flight. 

Under a given level of illumination, the pixel-to-pixel gain variation in the detector 
response generates artificial spatial fluctuations. The responsivity (flat-field) correction is 
critical to measuring spatial fluctuations in the DGL brightness. We measure the LRS flat- 
field response in the laboratory using the same set-up as described Subsection 13.21 

The measurement is made with two different integrating spheres, one with a 4 inch exit 
port and the other with an 8 inch port. The spheres are filled with light from three different 
light sources; the SIRCUS laser, the SCL and a quartz-halogen lamp with a solar-like filter 
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OZemcov et al. 2013). The variety of source-sphere combinations enables us to check the 
systematic uncertainty of the flat-held correction. The illumination patte rn of the output 
from these integrating spheres is uniform to better than 1% (Korngut et al. 2013), which is 
good tracer for the hat-held correction. 


Figure [6] shows the hat-held correction map. We calculated a hat hel, / = I\ >me &n/I\, P ixei, 
where /a, pixel indicates the brightness of the light sources detected by a pixel and /a, mean 
indicates the mean brightness of the light sources, for all pixels. Typically, the hat-held 
correction of / = 0.97-1.03 % applied to the image with an accuracy ~ 3 %. The accuracy 
of the hat-held is estimated from the variation between measurements with different set ups. 


4. Data Reduction 


4.1. Image from raw data 


The hrst step of the data reduction is to make an image from individual time-ordered 
array reads. For a charge integrating detector, the time derivative of the charge is propor¬ 
tional to the optical power . We ht a slope and offset to the raw data with a least squares 


method (Garnett & Forrest 19931 ). to derive the best-ht slope in photocurrent units (e- s 1 ). 


A spectral image of the NEP held from the second flight is shown in Figure [L] as an example. 
The hve vertical lines correspond to the images of hve slits, which are dispersed along the 
horizontal direction by the prism. Bright stars seen as bright horizontal lines are masked to 
obtain the diffuse signals only. The dark regions located at the top, bottom and left side are 
regions masked from infrared light. 


4.2. Removal of Bright Stars 

It is necessary to remove bright point sources from the image to isolate the diffuse 
components. We hrst average the photocurrent of each slit along the horizontal direction, 
then clip the pixels containing stars determined by the criterion that the band-averaged 
photocurrent is larger than the mean band-averaged photocurrent of all pixels by 2a, where 
a is the standard deviation of the photocurrent. We iterate this clipping procedure until the 
ratio of the number of rejected pixels to remaining pixels is less than 0.1 % of the total. An 
example masked image is shown in the right panel of Figure [0 All of the stars brighter than 
13th magnitude are removed by this procedure, and 95 % of 13th magnitude stars are also 
removed. 
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4.3. Dark Current 

In the absence of incident photons, the detector produces a small positive signal called 
“dark current”. In order to measure the absolute spectrum of the astrophysical sky, an 
accurate subtraction of the dark current is required. 

We estimate the dark current from the masked region of the array for each observation. 
The dark current is typically ~le- s -1 , corresponding to r^j 20 nW m 2 sr 1 at 1.25 /im. The 
dark current of the masked region is slightly different from that of the FOV region, based on 
shutter-closed data. This difference makes a ~ 0.03 e- s^ 1 systematic offset corresponding to 
0.7 nW m -2 sr -1 . As the DGL determination uses a fit to a spatial template, residual dark 
current adds a small amount of noise but does not bias the result as described in Section [7] 


5. Data Analysis 


The processed data includes astronomical emission from ZL, DGL, ISL and EBL, i.e., 


/ s k y — /zl + /dgl + hsh + I] 


EBL- 


( 1 ) 


We separate the DGL component using its spatial distribution as traced by 100 /an emission 
with structure on spatial scales smaller than a degree as shown in Figure [2] and |3] ZL 


is known to be spatially uniform on spatial scales smaller than a degree (1 Abraham et al. 


r 


19971 ; IK ash 1 in sky et al.l 120051: iPvo et al.l 120121 ) . In this analysis, we assume that there is no 
correlation between the EBL and the 100 /mi brightness. 

Although the small-scale fluctuations of the ZL brightness are negligible, the large- 
scale distribution may affect the correlation between the 100 /an brightness and the DGL 
brightness. To account for contamination from the ZL, LRS spectral images are constructed 
by aligning the DIRBE-based model predictions for large scale ZL structure with the location 
of the LRS slits. We use two different mo dels to check the systematic effect of the ZL 
subtraction ( Kelsall et ah 1998 ; Wright 1998 ). 


We first obtain a fiducial ZL spectrum by differencing the LRS data between holds, i.e, 
(/sky,i - hsL,i) ~ (/sky j - hsh,j)- The difference is approximately I Z L,i ~ Izlj when the DGL 
brightness is similar in held i and j, where Izl,i indicates the ZL spectrum, and I s ky,i indicates 
the sky spectrum. The ISL spectrum, /isl,u is estimated by Monte-Carlo simulations of 
the star distribution in the FOV using the 2MASS catalogue flSkrutskie et al.l [20061) and a 


population synthesis code for simulating Galactic star counts (IGirard i e t a l 


20051 


'), t aking 


into account the limiting magnitude and effective slit area of the LRS flMatsuura et ah 2014, 
in prep). We calculate the difference of every combination of helds and take an average of the 
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The NEP field of the second flight. 


The ElatlO field of the second flight. 
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Fig. 7. — The correlation of near-infrared surface brightness at 1.25 /jrn and 100 /mi bright¬ 
ness. Each panel shows a different field. There is a statistically significant correlation between 
near-infrared surface brightness and 100 /ini brightness in only these four fields. The small 
circles present the sky brightness, J s k y , with a best-fit line black solid line. The l arge circles 
indicate the subtracted ZL, Izl, which is estimated from a ZL model (iKelsall et al.l 19981 ), 
fitted as a dot-dashed line. The red asterisks indicates I s k y — Izl, which consists of DGL, 
1SL, and EBL, and is fitted by the dashed line. 
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The SWIRE field of the fourth flight. 
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Fig. 8.— Because the contrast of the DGL brightness is too low in these fields, there is 
no statistically significant correlation between near-infrared surface brightness and 100 /irn 
brightness. See Figure [7] caption. 
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Fig. 9.— The correlation of near-infrared and 100 //m brightness in the ElatlO field observed 
in the second flight. Each panel shows a different observed wavelength. See Figure [7j caption. 
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Fig. 10. 


Correlation in the NEP field observed in the fourth flight. See Figure [7] caption. 
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differences to make the ZL spectrum template. Because the difference of every combination 
is consistent with each other within the errors , we use this ZL spectrum template for every 
field. The a bsolute brightness of ZL is estim ated by using ZL models for the DIRBE data 
at 1.25 pm ( Kelsall et al. 1998 : Wright 1998). 


We separate the DGL component from J s k y — /zl using a linear correlation analysis. In 
the optically thin limit, the DGL brightness can be approximated to be: 

I DGL ~ 'YnirO'nirNIisrf, (2) 

where 7 nir and <Jnir present albedo and scattering cross sections of interstellar dust in the 
near-infrared respectively, Iisrf is the brightness of the ISRF, and N is column density of 
interstellar dust. In the optically thin limit, the far-infrared brightness is given by: 


I fir ~ (1 — , ') / fir)cifirNB(T), 


(3) 


where the B(T) is Planck function. We assume the temperature of interstellar dust is uni¬ 
form, so that the DGL brightness can be written as a function of the far-infrared brightness; 


l DGL 


Iniro nirIisrf 


l fir- 


(4) 


(1 — 'Jfir)^firB(T) 

Equation [4] indicates that the DGL brightness correlates linearly with the far-infrared bright 
ness in optically thin regions. 


Because thermal emission from interstellar dust dominates the diffuse sky brightness 
at 100 pm, it is easier to separate interstellar radiation from other diffuse sources in the 
far-infrared than in the near-infrared. As a result, we use a far-infrared intensity map 
measured at 100 pm ( Schlegel et al. 19981 ). The SFD 100 pm map is based on the all¬ 
sky survey combination of IRAS (IR Astronomical Satellite) and DIRBE/COBE (Diffuse 
Infrared Background Experiment / Cosmic Background Explorer). The SFD 100 pm map 
has the accurate calibration of DIRBE/COBE and the ~ 6' resolution of IRAS . 


Thus the linear correlation between our data and 100 pm brightness is written as: 

/sky AjL o(A) T 6(A)/i00/zm- (5) 

The slope 6(A) gives the conversion factor from 100 pm brightness to DGL brightness. The 
offset a(A) accounts for the ISL and EBL contributions. 


6. Results 

Figures [7] and [8] show the correlation between / s k y —/zl with the brightness of the 100 pm 
map in all fields at 1.25 pm. There is a statistically significant correlation in only four fields 
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with large spatial contrast in 100 /mi brightness, specifically the NEP held of the second and 
fourth flight, the DGL held, and the ElatlO held of the second flight. The Peason correlation 
coefficient, 


E -^data / T 

k =i 


r = 


100/im,k 


■ 100/im 


) (-^sky ,k 


sky 


'ZL ,fc 


I zl) 


E Nlata (T T \2. /V^Nlata IT T T 

k = i \-L100/j,m,k J-lOOfim) A/ / V-*sky,fc -l sky ^ZL ,k 


I zl) 1 


( 6 ) 


where A^ata indicates the number of data and / indicates the average of the data, is larger 
than 0.5 at A < 1.35pm only in these four helds. The Peason correlation coefficients of other 
helds are less than 0.5, indicating that the contrast of the DGL brightness is too low to 
measure DGL in other helds. Thus we use the results of the four statistically significant 
helds below. Figures [9] and [TO] show the correlation of / sky — J ZL with 100 pm brightness 
for the NEP held and the ElatlO held of the second flight at different wavelengths. We also 
present Izl, estimated from the ZL model and I sky in Figures 17118l 191 and flOl To improve S/N, 
array pixels are binned along the slit length direction and averaged into broad AA = 200 
nm wavelength bins. The error bars present 1 a statistical uncertainty estimated from the 
variance across all pixels in the binned region, added in quadrature with the ~3 % hat-held 
error. 


As shown in Figure |9] the scatter of J sky is traced by Izl, and the scatter of J sky — Izl is 
smaller than that of / sky . This result indicats that a large-scale ZL gradient is present, and 
that this gradient is explained by the ZL model. 

In Figures [7] O fTUl and l9l we ht the linear function from Equation [5] to our data, and 
Figure [III shows the slope b( A) as a function of wavelength. The results of the four helds are 
shown in Figure fill with the results of various optical and infrared measurements. 

The DGL spectrum of the NEP held observed in the second flight is consistent with 
that of the fourth flight within our errors. We obtain the same result from observations with 
different calibration measurements, which indicats that our measurements and calibrations 
are valid. The DGL spectrum of the ElatlO held is systematically low compared to the NEP 
held. 


Figure [13] shows the mean DGL spectrum of the four helds. DGL measurements at 
optical wave lengths indicate that the DGL brightness varies from cloud to cloud by a factor of 
3-4. Ienaka et ah (120131 ) argued that variation of optical depth and/or the forward scattering 
characteristics of dust grains might be the origin of this variation. Because we measured 
the DGL brightness in diffuse regions, away from dense molecular clouds and star forming 
regions, we only present optical measurements in similarly diffuse regions for comparison. 
However, because the DGL spectrum of BD12 had a large uncertainty in hux calibration, 
we do not present its spectrum. The Pioneer data (Matsu ok a et ah 2011]) are presented at 
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optical wavelengths, since the data covered one fourth of sky at galactic latitude higher than 
35°. As shown in Figure [13l our result is con sistent with interpolated levels from the Pioneer 
data and the AKARI data (Tsumura et al. 2013a). 


7. Systematic Uncertainty 

We estimate the possible systematic error from residual airglow emission, residual faint 
stars, and uncertainties from the ZL subtraction. 


7.1. Wavelength Calibration 

We estimate the systematic error from the wavelength calibration by shifting the wave¬ 
length map by ± 1 nm, which corresponds to the systematic error seen in multiple wavelength 
calibration measurements. We derive the DGL spectra using this shifted wavelength map 
and calculate the differences of the DGL spectra. As shown in Figure[[4l this mean difference 
on the mean DGL spectrum is negligible. 


7.2. Flat-Field Correction 

We bound the systematic error of the flat-field correction based on the difference be¬ 
tween the flat-field measurements using different integrating spheres. We use two integrating 
spheres with 4” and 8” exit port diameters as described in Subsection 13.31 We apply each 
flat-field and derive the DGL spectra in four fields. We then calculate the difference between 
the DGL spectra of each flat-field and calculate mean of this difference. We determine the 
mean difference as the systematic error on the mean DGL spectrum of the flat-field correc¬ 
tion, as shown in Figure HU The flat-field systematic error is < 6 %, which is acceptably 
small. 


7.3. Dark Current Subtraction 

To check the systematic error from the dark current subtraction, we use the shutter 
closed data acquired during the flights. We subtract the dark current from the shutter 
closed data using the same method for the field data described in 14.31 and then correlate 
with 100 /irn brightness. There is no correlation detected between the shutter closed data and 
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Fig. 11.— The DGL spectrum of all four fields. The intensity is normalized at I\ oo /4m = 1 
MJysr^ 1 . The red square indicates the DGL spectrum of the NEP field of the second flight, 
the cyan small circles the ElatlO held, the green diamonds the DGL held, and the blue 
triangles the NEP held of fourth flig ht. The error bars sh o w 1 a standard error. The opt i¬ 
cal measurements a re als o presented. lLaureiis et al.l (1198711 : iGuhathakurta fe Tvsonl (119891 1: 


Zag urv et al. (119991 1 and Paley et al. ( 1991 ) measure the DGL brightness correlating with 
the optical brightn ess w ith the ori gi nal IRAS J.00 /mi map. Recent optical measurements 

(12013111 correlate the optical bright- 


Ienaka et al. 


(jWitt et al.l (120081 1: iMatsuoka et al.l (1201111: _ 

ness with the SFD 100 fim map. Only Brandt fe Drama ( 2012 1 measures the spectrum of 
DGL at the optical wavelengths. Tsumura et ah (l2013bl ) measures the spectrum of DGL at 
1.8 ~ 5.0 n m and detects a polycyclic aromatic hydrocarbon feature. 
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wavelength [//m] 


Fig. 12.— Comparison between the mean DGL spectrum of all four fields, given by the solid 
line, and the ISL spectrum, given by the dashed line for stars fainter than 13th magnitude. 
The error bars give 1 a statistical error. 
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wavelength [pm] 


Fig. 13.— Comparison between the mean DGL spectrum and theoretical DGL models 
described in Table [2] The 100 /jrn intensity is normalized to that of 1 MJy sr' 1 . The 
red asterisks indicate our mean DGL spectru m given in Table [2l We only show the DGL 
brightness measured from diffuse sky regions flMatsuoka et al.l 1201 ll; iTsumura et al.l l2013bh 
to compare with our results. 
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Wavelength [//in] 


Fig. 14.— The systematic errors on the mean DGL spectrum. The intensity is normalized at 
hoo/im = 1 MJysr -1 . The instrumental systematic errors (purple line) associated with surface 
brightness calibration error (solid line), dark current subtraction (long dashed dotted line), 
flat-field error (dashed line), and wavelength calibration error (dotted line) are indicated. 
The systematic error from the astronomical foreground (cyan line) associated with airglow 
contamination (dotted line), ZL subtraction (dashed line), and residual faint stars (dotted 
line) are presented. The orange solid line indicates the total instrument and astronomical 
systematic error that is the quadrature sum of the systematic errors. The mean DGL spec¬ 
trum is presented as the red asterisks with la statistical error. An orange shaded region 
shows the total systematic error band about the mean DGL spectrum. 
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the 100 n m brightness; the Pearson correlation coefficient is less than 0.3 in all wavelengths 
at all helds. The mean slope of this correlation of four fields is determined to be the worst- 
case systematic error of the dark current subtraction and is shown in Figure [TT] as a function 
of wavelength. 


7.4. Airglow 

To estimate the contamination from airglow, we extracte the spatial structure of the 
airglow emission by differencing the first-half and the second-half integration of each held. We 
attribute a time and altitude dependence to the airglow emission, so the observed brightness 
is written as 

-fob s (t, h ) = J sky + / air (t, h) (7) 

where I air indicates the brightness of airglow emission, t indicates time from the launch, and 
h is altitude. We make an image of the spatial structure of airglow by differencing the first- 

half and the second-half integration, -f 0 bs(f first) hfirst) -fobs {t second) ^second) -fair (f first) f firsl. ) 
-fair(fsecond) ^second)) to cancel the astronomical component. 

To estimate the systematic contribution of the airglow emission to the DGL spectrum, 
we correlate this airglow image with the 100 /mi brightness map. The Pearson correlation 
coefficient is less than 0.5 in all wavelengths at all holds. This result indicates that there is 
no correlation between the airglow emission and 100 /mi brightness. From these results, we 
conclude that the airglow emission does not have a systematic effect on the DGL measure¬ 
ment. We estimate the maximum possible systematic error of airglow on the mean DGL 
spectrum from the mean slope of this correlation between airglow image with the 100 /mi 
brightness map of four helds, as shown in Figure [HJ 


7.5. Zodiacal Light Subtraction 


The systematic uncertainty of the ZL subtraction is quantified by calculating the dif¬ 
ferenc e between the DG L spectrum de rived using two different ZL models for the COBE 
data ( Kelsall et ah 1998 ; Wright 1998). In the ElatlO held, because the direction of the 
large-scale ZL gradient is similar to the spatial structure of DGL as shown in Figure [9j the 
systematic uncertainty of the ZL subtraction is 10 ~ 30 % of the DGL brightness, depending 
on the wavelength. In the other four helds, the systematic uncertainty of the ZL subtraction 
is ~1 % of the DGL brightness. The systematic error of the ZL subtraction on the mean 
DGL spectrum is presented in Figure dH 
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7.6. Residual Faint Stars 

Residual faint stars may spatially correlate with 100 /im brightness and make a system¬ 
atic effect. To rule this out, we correlate the integrated brightness of stars fainter than 13th 
magnitude from 2MASS at H- and J-bands with the 100 /im brightness map. Although the 
contribution of faint stars is brightest in the NEP fields, the Pearson correlation coefficient 
of the NEP fields is ~ 0.3 at H and J-bands. This indicates that there is no correlation 
between the faint stars with the 100 pm brightness. There are also no correlations in the 
DGL and Elat 10 fields. We determine the worst-case systematic error on the mean DGL 
spectrum as the mean slope of this correlation between faint stars and 100 pm brightness 
map of four fields, as shown in Figure ITT 


7.7. Contribution of Extra-Galactic Background Light 


According to S chleg el et al. (119981 ). galaxies with fluxes brighter 1.2 Jy were removed 
when they construct the SFD 100 pm map. Gala xies fainter than this limit have a peak 
flux density of ~ 0.12 MJy sr -1 at 100 pm. ISchlegel et al.l (ll998l) . however, mentioned 
that the contamination from extragalactic objects was very nearly uniformly distributed. 
Matsuura et al.l (20111) measured the far-infrared sky by AKARI at 90 pm and estimated 
that the spatial fluctuation due to galaxies was < 0.01 MJy sr -1 at 1° scale before removing 
bright galaxies. Thus galaxies contributed < 1 % of the DGL measurement, since the spatial 
fluctuation of the SFD 100 pm map is > 1 MJy sr -1 in our fields. We show upper limit to 
the systematic error from residual galaxies on the mean DGL spectrum in Figure dH which 
is negligible for our DGL measurements. 


7.8. Non-linear correlation between DGL and 100 pm brightness 


If the dust scattering is_optica 


lyjdiick, near-infrared DGL will not correlate linearly with 


100 pm brightness. Ien ak a et al. (120 13 ) presented the linear correlation breaks at > 6 MJy 
in the optical wavelengths. These results indicated that the linear correlation b reaks 


sr 


appe ar at extinction A(A) > 0.5, where A(A) is extinction at wavelength A (ISchle gel et_ al. 


1998 ). Since the 100 pm brightnesses in our observed fields is always < 6 MJy 1 correspond¬ 


ing to A(0.9 pm) = 0.4, we expect the correlation to be linear. 
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7.9. Systematic Error of the Mean DGL Spectrum 

We combine the systematic errors of the final mean DGL spectrum as shown in the 
orange band of Figure [14j The total systematic error is 30 % of the mean DGL spectrum 
at A < 1.35 pm, and > 30 % of the mean DGL spectrum at A > 1.35 pm. The total 
systematic error is the maximum case because there is no statistically significant correlation 
between dark current, airglow and residual faint stars with 100 pm brightness. 


8. Discussion 


As seen in Figure HU the slope 6(A) of the NEP and DGL fields are consistent. The 
ElatlO held is marginally lower in amplitude, although the errors are large, ft is possible that 
the large ZL correction in this held is partly responsible, or the result may be from spatial 
dependence of the DGL spectrum. Figure 4 of BD12 also shows the Galactic latitude and 
the Galactic longitude dependence of the slope 6(A) at optical wavelengths. If the difference 
in slope of 6(A) between helds is real, it could be attributed to spatial variations of dust 
scattering properties in the ISRF spectrum. However, it is difficult to claim significance 
from our data alone. To definitively detect spatial variations, higher S/N observations and 
wider sky coverage are required. 


Figure [12] shows the mean spectrum of stars fainter than 13 magnitude in the NEP 
held, obtained from the TRILEGAL code (j Girardi et all 120051 ). The mean DGL spectrum 
shows a bluer color than the ISL spectrum. Small grains are mainly responsible for DGL 
production, since scattering by large particles would not show the observed reddening. The 
LRS result is consistent with measurements of polarization of starlight by interstellar dust, 
which suggest the gra ins responsible for polarization of starlight have diameter a ~ 0.1 pm 
fjKim fe Martin] Il995l) . This result is also consis tent with observed stellar ex t inction from 
the i nfrared, 4 pm, to the ultraviolet, 0.1 pm (jWeingartner fe Draine 2001 : Zubko et al. 


20041 ). WD01 and ZDA04 constructed models to reproduce the observed stellar extinction, 
and argued that interstellar dust is composed of grains with 10~ 3 pm < a < 1 pm peak at 
a ~ 0.1 pm. 


To constrain the size distribution of dust particles from the LRS measurements, we 
compare the mean DGL spectrum with theoretical models from BD12 in Figure [T3l BD12 
considered an infinite plane-parallel Galaxy with a G aussian vertical distribution of dust from 
cr = 25pc flMalhotral 1 19951: iNakanishi fe Sofuel 120031). and a two-exponential distribution of 
stars with scale heights of 300pc and 1350pc (IBinnev fe Merrifieldl Il998l: iGilmore fe Reid 


19831) . BD12 estimated the stellar emission spectrum in two ways: (1) a model that repro- 
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duces the local ISRF of Mathis et al. 11 983 1 (hereafter MMP83), and (2) a stellar population 


synthesis model of Bruzual & Chariot (120031 ) (hereafter BC03) with solar metallicity and an 
exponential star formation history over 12Gyr. BC12 also used two typical dust models from 
WD01 and ZDA04. The dust of the WD01 model consisted of graphite, silicate, and PAH 
material, while the dust composition of the ZDA04 model consisted of bare graphite grains, 
bare silicate grains and PAHs. The size distribution of the ZDA04 model was shifted to 
smaller grains compared with WD01 model. The half-mass grain radius was a 0 .5 ~ 0.12 /un 
for both silicate and carbonaceous grains in the WD01 model, and dog ~ 0.06 /ini and 


0.07 /irn for carbonaceous grains and silicate grains respectively (ID raine 201 ll). Dust grains 


of size a > 0.2 fin i were absent in the ZDA04 model. We summarize these models in Ta¬ 
ble |2j Because the mo dels underestima t e b( A) by a factor of 2, the model is normalized at 
0.44 /iin to the data of Matsuoka et al. (20111). Iena ka et ah (120131 ) presented two possible 
explanations for this; a deficiency in UV photons in the ISRF, or that the dust grain albedo 
is different higher than what is assumed in the models. 


The combination of our result, the Pioneer data, and the shortest band of the AKARI 
data can be reproduced by both models of ZDA04. On the other hand, the models of WD01 
cannot reproduce the observed DGL spectrum. The models do not take into account PAH 
emission and thermal emission, so they do not reproduce the DGL spectrum measured by 
AKARI. Because the effect of different interstellar radiation fields is not significant at near- 
infrared wavelengths, our result implies that interstellar dust is dominated by small particles 
with Go .5 ~ 0.06 /irn, with a few large grains a > 0.2 /irn. 


Table 2. A summary of the interstellar dust scattering models of Brandt & Drain (2012), 
based on Weingartner & Draine (2001) (WD01), Zubko et al (2004) (ZDA04), Maths et al. 
(1983) (MMP83), and Bruzual & Chariot (2003) (BC03). 


Model Name 

Composition 

Interstellar Dust Size 

Interstellar Radiation Field 

WDCU-BC03 

WD01-MMP83 

ZDA04-BC03 

ZDA04-MMP83 

Graphite, silicate, and PAH 
Graphite, silicate, and PAH 
Bare graphite, bare silicate, and PAH 
Bare graphite, bare silicate, and PAH 

3 . 0 . 5 = 0.12 

3 - 0 . 5 — 0-12 

3 . 0 . 5 — 0.06~0.07 
3 - 0 . 5 — 0-06~0.07 

Stellar population synthesis 
Solar neighborhood 
Stellar population synthesis 
Solar neighborhood 
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9. Summary 

We measure the spectrum of diffuse galactic light in the near-infrared, helping to de¬ 
termine the properties of interstellar dust, particularly its size distribution. Since airglow 
emission is too bright to measure the DGL spectrum from the ground, the spectrum has 
never been observed at diffuse sky regions in the near-infrared. To derive the DGL spectrum 
as shown in Figure HU we correlate spectral images measured by the LRS with 100 /ini 
intensity. The measured DGL spectrum shows no resolved spectral features and is smoothly 
connected to the other DGL measurements in the optical and near-infrared wavelengths 
longer than 1.8 /irn, as shown in Figure ITol Rayleigh scattering of starlight by small grains 
largely explains the DGL spectrum. Our result implies that the size distribution is composed 
of small grains with a half-mass grain radius ao .5 ~ 0.06pm. 
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